QND Measurement of Large-Spin Ensembles by Dynamical Decoupling 
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Quantum non-demolition (QND) measurement of collective variables by off-resonant optical prob- 
ing has the ability to create entanglement and squeezing in atomic ensembles. Until now, this tech- 
nique has been applied to real or effective spin one-half systems. We show theoretically that the 
build-up of Raman coherence prevents the naive application of this technique to larger spin atoms, 
but that dynamical decoupling can be used to recover the ideal QND behavior. We experimentally 
demonstrate dynamical decoupling by using a two-polarization probing technique. The decoupled 
QND measurement achieves a sensitivity 5.7(6) dB better than the spin projection noise. 
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Quantum non-demolition measurement plays a cen- 
tral role in quantum networking and quantum metrol- 
ogy for its ability to simultaneously detect and generate 
non-classical quantum states. The original proposal by 
Braginsky [-] in the context of gravitational wave detec- 
tion has been generalized to the optical [2, ], atomic [ ] 
and nano-mechanical [ ] domains. In the atomic domain, 
QND by dispersive optical probing of spins or pseudo- 
spins has been demonstrated using ensembles of cold 
atoms on a clock transition [« , ], and with polarization 
variables [ , ], but thus far only with real or effective 
spin- 1/2 systems. 

QND measurement of larger spin systems offers a 
metrological advantage, e.g., in magnetometry [^^^], and 
may be essential for the detection of different quantum 
phases of degenerate atomic gases that intrinsically rely 
on large-spin systems [ - ]. Dispersive interactions 
with large-spin atoms are complicated by the presence 
of non-QND-type terms in the effective Hamiltonian de- 
scribing the interaction [ - ]. As we show, and con- 
trary to what has often been assumed [ -13, 17], these 
terms spoil the QND performance, even in the large- 
detuning limit. The non-QND terms introduce noise into 
the measured variable, or equivalent ly decoherence into 
the atomic state. The problem is serious for both large 
and small ensembles, so that naive application of disper- 
sive probing fails for several of the above-cited proposals. 

We approach this problem using the methods of dy- 
namical decoupling [ 3-zu], which allow us to effectively 
cancel the non-QND terms in the Hamiltonian while re- 
taining the QND term. To our knowledge, this is the first 
application of this method to quantum non-demolition 
measurements. Dynamical decoupling has been exten- 
sively applied in magnetic resonance [21, 22], used to 
suppress collisional decoherence in a thermal vapor [ ], 
to extend coherence times in solids [ ], in Rydberg 
atoms [ ], and with photon polarization [26]. Other 
approaches include application of a static perturbation 



[27, 28]. 

We consider an ensemble of spin-/ atoms interacting 
with a pulse of near-resonant polarized light. As de- 
scribed in references [14-16], the light and atoms interact 
by the effective Hamiltonian i^eff 



Ti^eff = GiSzJz + G2(5'xJx + *5'yJy) , 



(1) 



where r is the duration of the pulse and Gi^2 are coupling 
constants that depend on the atomic absorption cross sec- 
tion, the beam geometry, the detuning from resonance 
A, and the hyperfine structure of the atom [ ]. The 
atomic variables J (described below) are collective spin 
and alignment operators. The light is described by the 
Stokes operators S defined as Si = ^(^!^, aL)cri(a+, a_)-^, 
where the a^ are the Pauli matrices and d± are annihi- 
lation operators for the temporal mode of the pulse and 
circular plus/minus polarization. Bold subscripts, e.g., 
X, are used to label non-spatial directions for atomic and 
light variables. The Gi term describes a QND interac- 
tion, while the G2 describes a more complicated coupling. 
In the dispersive, i.e. far-detuned, regime, Gi and G2 
scale as A~^ and A~^, respectively. It has sometimes 
been assumed that the G2 terms can be neglected for 
sufficiently large A, leaving an approximate QND inter- 
action. As we show below, this scaling argument fails, 
and the G2 terms remain important. We note an impor- 
tant symmetry: i/^eff commutes with Sz-\- Jz, and is thus 
invariant under simultaneous rotation of J and S about 
the z axis. 

The atomic collective variables are Jk = ^^ ^ j^ 
where the superscript indicates the i-th. atom and jx = 
(/2-/2)/2, jy ^ (/./, + /,/.)/2, j, ^ /,/2 and j[x,y] ^ 

-^[jx, jy] = fzW^ - fl - 1/2). These obey commutation 
relations [jz,jx] = Uy, [jy,Jz] = Ux, [jx,jy] = U[x,y]. 
For / = 1/2, jx7 jy and j[x,y] vanish identically while 
for / = 1, j[x,y] = jz so that jx, jy, and j^ describe a 
pseudo-spin j. 



In the QND scenario, an initial coherent polarization 
state with (S) = (7Vi,/2,0,0) is passed through the en- 
semble and experiences a rotation due to the Gi term 
such that the component Sy (the 'meter' variable) indi- 
cates the value of J^ (the 'system' variable). We assume 
that Jx = Na/2. For a weak pulse, i.e., for (S) suffi- 
ciently small, we have the r-linear input-output relations 
^(out) ^ ^(in) _ ir[A^^^\Hefi]. Of specific interest are 
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(2) 
(3) 

(4) 



which describe the change in the system variable, its con- 
jugate, and the meter variable. In the case of / = 1/2, 
the G2 terms vanish identically and we have a pure QND 
measurement: information about Jz enters ^y and there 
is a back-action on Jy, but not on Jz. The input noise 
var(6'y ^) = 6'x/2 limits the performance of the mea- 
surement, and corresponds to a spin sensitivity of SJ^ = 
(2G^5'x)~^. For comparison, the projection noise of an x- 
polarized spin state is var(Jz) = Jx/2, so that projection 
noise sensitivity is achieved for S'x = (G^i Jx)~^ = *5'snr- 

This ideal QND regime does not occur naturally ex- 
cept for / = 1/2. In the interesting regime ^x ^ -S'snr, 
we find that G2Sy^Jy ^ Jy{G2/G\)/ Jy^ is independent of 
A, and cannot be neglected based on detuning. To get 
an order of magnitude, we note that for large detuning, 
Gi ^ (JqT /AA/S.^ G2 ^ GiAhfs/^ where ctq is the on- 
resonance scattering cross-section, A is the effective area 
of the beam, and F and Ahfs are the natural linewidth 
and hyperfine splitting, respectively, of the excited states. 
In terms of the on-resonance optical depth do = ctqA^a/^, 
we find G2/G\Jy^ ~ SAhps/^^oT. In a typical experiment 
with rubidium on the D2 line, Anrs/r ~ 30 and d^ ~ 50 
[29], so the contribution of this term is important. 

In contrast, the last term in Eq. (3) and (4), respec- 
tively, contribute variances (G^SyJ^) and 

which scale as A~^. We will henceforth drop these terms. 
The system variable Jz is coupled to a degree of free- 
dom, Jy, which is neither system nor meter in the QND 
measurement. This coupling introduces noise into the 
system variable, and decoherence into the state of the 
ensemble. To remove the decoherence associated with 
this coupling G2Sy^Jy^ we adopt the strategy of "bang- 
bang" dynamical decoupling [18-20]. In this method, a 
unitary [/^ and its inverse Ul are alternately and period- 
ically applied to the system p times during the evolution, 
so that the total evolution is [Ul(jH(t/2p)tJi,tjH(t/2p)Y 
where Unit) describes unitary evolution under H for a 
time t. With this evolution, those system variables that 
are unchanged by Ub continue to evolve under iJ, while 



r^2 q2 t2 
'-^2^z^y 



others are rapidly switched from one value to another, 
preventing coherent evolution. For large p, the system 
evolves under a modified Hamiltonian H' = PH^ where 
P projects onto the commutant (i.e., the set of operators 
which commute with) of {[/fo, f/J} [20]. 

To eliminate G2(5'xJx + 5'y Jy), while keeping GiSzJz 
we choose a t/5 which commutes with J^, but not with 
Jx or Jy, namely a tt rotation about Jz, Ub = exp['i7rJz]. 
This leaves Jz unchanged, but inverts Jx and Jy. By the 
symmetry of iJeff? Ihis is equivalent to inverting ^x and 
Sy^ which suggests a practical implementation: probe 
with pulses of alternating S'x, and define a 'meter' vari- 
able taking into account the inversion of ^y. 

We consider sequential interaction of the ensemble with 
a pair of pulses, with S^ ^ = —S^ = Nl/4p. We de- 
fine also the new 'meter' variable Sy ^ = Sy — Sy . 
We describe the atomic variables before, between, and 
after the two pulses with superscripts (in), (mid), (out), 
respectively. We apply Equations (2-4) to find: 
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and 



plus terms in GiG25'x5'z Jx, ^2^x^[x,y] and GiG2S^Jy 
which become negligible in the limit of large p. The ideal 
QND form is recovered by the dynamical decoupling. 

The presence of the G2 term can be detected by noise 
scaling properties. While in the ideal QND of Equations 
(8), (9) the variance of the system variable is ex Jx giving 
a variance for the meter variable linear in Jx, for the 
imperfect QND of Equations (2) to (4) this is not the 
case: from Equation (6), we see that Jy acquires a back- 
action variance oc J^, which then is fed into the system 
variable by the G2 term. This additional J^ noise is also 
reflected in the meter variable, and provides a measurable 
indication of G2 . 

We use the two-polarization decoupling technique to 
perform QND measurement on an ensemble of '^ 10^ 
laser cooled ^^Rb atoms in the F = 1 ground state. In 
the atomic ensemble system, described in detail in ref- 
erence [ ], /is pulses interact with an elongated atomic 
cloud and are detected by a shot-noise-limited polarime- 
ter. The experiment achieves projection noise limited 
sensitivity, as calibrated against a thermal spin state [ ] . 

The experimental sequence is shown schematically in 
Fig. 1. In each measurement cycle the atom number 
Na is first measured by a dispersive atom-number mea- 
surement (DANM) [ ]. A Jx-polarized coherent spin 
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Figure 1: (color online) Experimental sequence for projection 
noise measurement. The CSS is prepared once and its magni- 
tude ( Jx ) is measured. This serves as a measure of the spin 
polarization prior to the QND probing. We prepare the CSS 
a second time and assume it has the same spin polarization 
as in the first preparation. The state is probed with a train 
of pulses of alternating polarization. Measuring the spin po- 
larization after the QND measurement tells us the amount 
of depolarization introduced in the QND probing. The QND 
probing scatters a non- negligible fraction of atoms into F — 2, 
which are removed from the trap with resonant light in order 
to reduce the number of atoms in the trap. The whole cycle 
is repeated 10 times during one trap loading. 



state (CSS) is then prepared and probed with pulses of 
alternating polarization to find the QND signal ^y = 



E, 



g(out)/ 



( — 1)*+^. Immediately after, ( Jx) is measured 



to quantify depolarization of the sample and any atoms 
having made transitions to the F = 2 manifold are re- 
moved from the trap, reducing Na for the next cycle and 
allowing a range of Na to be probed on a single load- 
ing. This sequence of state preparation and probing is 
repeated ten times for each loading of the trap. The trap 
is loaded 350 times to acquire statistics. 

The optical dipole trap, formed by a weakly- focused 
(52 /im) beam of a YbiYAG laser at 1030 nm with 6W 
of optical power, is loaded from a conventional two 
stage magneto-optical trap (MOT) during 4 s. Sub- 
Doppler cooling produces atom temperatures down to 
25 /iK as measured in the dipole trap [29]. In the 
DANM, we prepare a Jx-polarized CSS, i.e., all atoms 
in a coherent superposition of hyperfine states |t / i) = 
|F = 1, mp = ±1), by optically pumping with vertically- 
polarized light tuned to the transition F = 1 ^ F' = 1^ 
while also applying repumping on the F = 2 ^ F' = 2 
transition and a weak magnetic field along x to prevent 
spin precession. The atoms arrive to this dark state after 
scattering fewer than two photons on average. To mea- 
sure ( Jx), we send ten circularly-polarized probe pulses, 
i.e., with (5'z ) = Nl/2, tuned 190 MHz to the red of the 
transition F = 1 ^ F' = {). Each pulse, of 1 jis dura- 
tion, contains 2.6 x 10^ photons and produces a signal 
(6'y) (X G2 (^Sz) ( Jx). The coherent state for the QND 
measurement is prepared in the same way, but in zero 
magnetic field. 

To measure Jz, i.e., one half the population difference 
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Figure 2: (color online) Variance of polarimeter signal as a 
function of atom number, comparing naive probing, i.e., a sin- 
gle input polarization, to "bang-bang" dynamically-decoupled 
probing of different orders p. Grey curves indicate simula- 
tion results for: naive probing (solid), and decoupled probing 
with ip — \ (widely dashed), p — 2 (dashed), and p = 5 
(dotted). The black solid line shows the expected projec- 
tion noise for p ^ 00, or the ideal QND interaction G2 — 0. 
All curves are calculated using the independently measured 
interaction strength G\ — 1.27(5) x 10"^ and have no free 
parameters. Red squares are measured data using dynamical 
decoupling with p = 5. Blue circles are measured data with 
naive probing. Technical noise from laboratory fields domi- 
nates the naive probing results, and pushes them above the 
theoretical curve, while technical noise is suppressed in the 
dynamically-decoupled probing. 



between \X) and ||), we send probe pulses of either ver- 
tical 5x = Til, 12 or horizontal 5x = —71^/2 polarization 
through atomic sample and record their polarization ro- 



tation as Sy ^ . The number of individual probe pulses 
is 2p and the total number of probe photons Nj, = 2pnL . 

In Fig. 2 we plot the measured noise versus atom num- 
ber, which confirms the linear scaling characteristic of 
the QND measurement. The black squares indicate the 
variance var(S'y) normalized to the optical polarization 
noise, measured in the absence of atoms. Independent 
measurements confirm the polarimetry is shot-noise lim- 
ited in this regime. The black solid line is the expected 
project ion noise scaling 4 var(6'y)/A^L = l+GiA/'Lvar(Jz), 
calculated from the independently measured interaction 
strength Gi and number of probe photons Nl = 8 x 10^. 
The QND measurement achieves projection- noise limited 
sensitivity, i.e., the measurement noise is 5.7(6) dB below 
the projection noise. 

Also shown are results of covariance matrix calcula- 
tions, following the techniques of reference [^v.], includ- 
ing loss and photon scattering. The scenarios considered 
include the naive QND measurement, i.e., with a sin- 
gle polarization, and the "bang-bang" or two-polarization 



QND measurement, with p — 1,2,5. These show a rapid 
decrease in the quadratic component with increasing p. 
This confirms the removal of G2 due to the dynamical 
decoupling. Also included in these simulations is the 
term SyJy which introduces noise into Jz proportional 
to G2var(S'y) ( Jx) . For our experimental parameters 
this term leads to an increase of var( Jz) of less then 2 % 
and as noted above could be reduced with increased de- 
tuning. 

The dynamical decoupling also suppresses technical 
noise which would otherwise enter into Jz through the 
interaction G2{S^J^ + ^y-^y)- An imperfect preparation 
of the atomic and and/or light state, e.g., ( Jy) / or 
(^y) 7^ 0, would otherwise be transferred into Jz. 

Using dynamical decoupling techniques, we have 
demonstrated optical quantum non-demolition measure- 
ment of a large-spin system. We first identify an often- 
overlooked impediment to this goal: the tensorial polariz- 
ability causes decoherence of the measured variable, and 
prevents (naive) QND measurement of small ensembles. 
We then identify an appropriate dynamical decoupling 
strategy to cancel the tensorial components of the effec- 
tive Hamiltonian, and implement the strategy with an 
ensemble of ~ 10^ cold ^^Rb atoms and two-polarization 
probing. The dynamically-decoupled QND measurement 
achieves a sensitivity 5.7(6) dB better than the projec- 
tion noise level. The technique will enable the use of 
large-spin ensembles in quantum metrology and quan- 
tum networking, and permit the QND measurement of 
exotic phases of large-spin condensed atomic gases. 
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